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We study a model of nonequilibrium quantum transport of particles and energy in a 
many-body system connected to mesoscopic Fermi reservoirs (the so-called meso-reservoirs) . 
We discuss the conservation laws of particles and energy within our setup as well as the 
transport properties of quasi-periodic and disordered chains. 

§1. Introduction 

Understanding the macroscopic transport from a microscopic point of view is 
a central topic of statistical physics. Especially, the development of nanoscale de- 
vices reveals unconventional transport. Thus, the study of transport in mesoscopic 
systems has significant value both in fundamental theory and for applications in 
future technologies. Since mesoscopic systems are strongly coupled with the envi- 
ronment, their understanding in nonequilibrium regimes requires knowledge of the 
global features of the total system including the environment (reservoirs). Recently, 
we proposed a model that comprises mesoscopic reservoirs with a finite number of 
degrees of freedom^ (the so-called meso-reservoirs), and study the transport prop- 
erties of periodic chains. In particular, the parameter dependence of the transport 
properties as well as the Onsager reciprocity relation were studied. 

In this paper, we study periodic chains in more detail, as well as quasi-periodic 
and disordered chains. We will also discuss the conservation law of energy and 
particles within our setup. 

§2. Model 

A key idea of our model is to enforce the finite reservoirs to equilibrium (or almost 
equlibirum) state using the Lindblad dissipatorP'® In particular, if the term in the 
Liouvillean evolution containing the Lindblad dissipator is small, we can interpret 
that the 'traced-out-infinite-reservoirs' {super-reservoirs) force our finite reservoirs 
{meso-reservoirs) to equilibrium. That is to say, our density matrix follows the 
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Lindblad equation of the following form: 
dp 



- i [H,p] + D(p) (2-1) 

D (P) = Yl ( 2L k^mpLl ^ m -{Ll a m L k ^ m ,p} 

k,a,m 

H = Hs + H L + H R + V 

n—1 n 

H s = -J2 • ' + ( h - c -)) + E U J c fe 

3=1 3=1 
K 

H a = y^ y e%al a a ka , a = L,R 

k=l 
K 

V = Y, (^4l C 1 + V Hn C n) + (h.C.) 
k=l 

Lk,a,l = sj A,o,lflfca) ijfe,a,2 = V^\a^2 a la 

A,a,i = 7fc(l-^a(e fc )), r feja>2 = j%F a (e k ) , 



where Cj is the annihilation operator of system fermions, a k: a is that of reservoir 
fermions with wave number k, tj is the nearest neighbor hopping, Uj is the on- 
site potential, v? is the coupling between the system and reservoirs, and F a {e) = 
( e A*(e-/A*) _|_ are F erm i distributions, with inverse temperatures j3 a and chemical 
potentials p a , while [-, •] and {•, •} denote the commutator and anti-commutator, re- 
spectively. The parameter 7^ determines the strength of the coupling to the super- 
reservoirs and needs to be fine-tuned in order to ensure the applicability of the 
modelp} We stress that our model does not rely on the usual weak-coupling as- 
sumption needed for the physical derivation of the Lindblad master equation,® thus 
7^ do not need to be small parameters. 

§3. Conservation laws 

3.1. Particle conservation 

In this section, we discuss conservation of particle number and energy. We focus 
first on the particle current. Let nj = S-Cj and = a\ a a ka (a = L,R), then, with 
time derivatives that can be casted as 

Jf-x ~ Jf, (1 < 3 < n - 1) 

*X^Ki C l - C Wl) = ^3k 
k k 

ilji'Yj ■ : - c ]+i c i)' (1 — J — n — 1) 

-«X^(4il C n ~ c i a kR) = ^2 3k 
k k 



dnj 
~dt ~ 



jP 

■J , , 
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d 4 -jf 

^=jf + D(n«) , (3-1) 

where the dissipative parts D{n k ) are given by 

£>("*) = -2r k>a>1 n a k + 2r fe)Qi2 (l - r%) 

= -2^n a k +2^F a (k)l . (3-2) 

By taking the expectation value in the nonequilibrium steady state (NESS) (•} = 
tr(-p(i — > oo)), we obtain 



= < J °> = (Jl) = - = (Jn-l) = (Jn) = ) = J P 

k k 

(j£) = (D(n L k )), Of ) = -{D(n£ )) . (3-3) 



Thus, the particle current from the left meso-reservoir to the system ^2 k (ju) is equal 
to the current from the system to the right meso-reservoir ^2 k (j k ), and there is no 
particle loss due to the existence of the Lindblad dissipators. Moreover, we have 

J P = ^2(D(n L k )) = -J2(D(nf)), (3-4) 

k k 

and it follows 

J p = 2^ 7fc L {(n L k ) - F L (k)} = -2£t£ ? {(nf ) - F R (k)} . (3-5) 

k k 

This expression can be understood as a generalization of the Landauer formula for 
the current between the meso-reservoirs and the super-reservoirs. Moreover, if j k 
is independent of k, the total differences of charge densities from their equilibrium 
state values are proportional to (1/7°). 

3.2. Energy conservation 

We now turn our attention to the conservation of energy. Casting the local 
energy density as 

Hj = -{t jC )c j+ i + t jC ) +1 c^ + U jC ] Cj , (1 < j < n) 
c n +i = , (3-6) 

their time derivatives are 
dH 



-^ = jf_ 1 -jf, (1<j<»-2) 
Jq = UiJ p + i^v L ti(c\a k - a\c 2 ) 
jf = itjtj. \{r ■(■;.■> ~ 4+2 c i) + u j+i J fi (1 < 3 < n ~ 2 ) 
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^ ~ J n-2 J n-1 J n 

jE = rr jP 
J n-1 — u n<J n -l 



J n = UJ n = iU Yl V k( C i a kR ~ a kR C n) 



k 

T E T E 



dV L 
dt 
dV R 



k 

tE , T E 



-J£ + J^ V + D(V L ) 



dt 



J^n = -iJ2 e kVk(b\c n - c n b k ]) = ^ 4j{ 
dH L 



dt 
dH R 



-JE^ V + D(H L ) 

J^ R + D{H R ) , (3-7) 



dt 

where the dissipative terms D(V a ) and D{H a ) are given by 

k 

D(V R ) = -Y,^kVk{A R C n + 4flkR) 
k 

D{H a ) = 2Y J l%tk(F L {k)l-nt) . (3-8) 

k 

By taking average with respect to the NESS we obtain 

(Jf > = (Jl) = ■■■ = (Jn-2) = (Jn-l) + (Jn) = ^n) + (%) 

»J^«fc*l(4 a * - 4 C 2) = tn-lVk(a\ R C n -l ~ 4_l a fcfl) 

k k 

(D(H L )) + (D(V L )) = —(D(H R )) - (D(V R )) . (3-9) 

Although, the total particle current at the two meso-reservoirs is conserved, i.e., 
J2kUk ) = Y!k(jk)> ^ ne particle current distributions are not the same, i.e., 

- 2^{(n L k ) - F L (k)} = <#) + (jf) = 2 7 *{<n£) - F R (k)} . (3-10) 

Thus, the ingoing energy current at the left hand sice is not equal to the outgoing 
energy current at the right hand side, i.e., Ylk e k(^k) ^ ^2k e ktik)- 
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In our case, 

(D(H L )) = J>(j fc L ) = 2J2^ IMcW) (3-11) 
k k 

and 

(D(V L )) = -2j2lk4 McW) (3-12) 
k 

are not conseved separately, and it induces the difference of the energy current from 
left meso-reservoir to system and that from system to right meso-reservoir. Namely, 
the following four statements are deeply connected: 

(1) Different amount of energy is dissipated at left and right (D(Vl)) ^ — (D(Vr)). 

(2) The energy current from left meso-reservoir and that to right meso-reservoir 
are different (J^+ s ) + {Js^r)- 

(3) The particle current distribution of two reservoirs are different (jf) ^ (Jk)- 

(4) The particle distribution of two reservoirs are different (n^) ^ (n^). 

§4. Numerical results 

As discussed inP expectation values with respect to the NESS are easily ob- 
tained by solving 2n x 2n dimensional Sylvester equation. In this section, we discuss 
the transport properties of periodic, Fibonacci and disordered chains at NESS. 

We study monoatomic (tj = t), diatomic (£2.7-1 = tA, = ts), Fibonacci, 
and disordered chains. The Fibonacci chain is constructed by first taking n = 3 
and setting t\ = tA, £2 = £b constituting the first generation, and then inductively 
replacing tA by tAtB and ts by tA- For instance, the second generation yields 
n = 4 and t\ = tA, ti = £g, £3 = tA', the third generation yields n = 6 and 
t\ = tA, t2 = ts, ta = tA, ti = tA, £5 = ts, and so on. Alternatively, one can 
construct n-th generation by concatenating (n — 2)-th generation after (n — l)-th 
generation. 

For disordered chains, we take tj from a uniform distribution in [t — S, t + 6]. 

To make things simple, we have set Uj = U, v% = vj} = v and 7^ = 7^ = 7 
for numerical results. As discussed should be smaller than 7, and we have 

set 7 = 0.1, v = 0.03, ei = —20, ex = 20, K = 200, and hl = —fJ-R = V, unless 
specified differently 

We have checked that reservoirs are close enough to equilibrium and satisfy 
conservation law (|3-5|) . We show the j dependence of occupation density (rij) and 
n dependence of the particle current J p in Fig. [1] (monoatomic t = 3 and diatomic 
chain tA = 3, ts = 6), Fig. [2] (Fibonacci chain tA = 3, ts = 6), and Fig. [3] (disordered 
chain t = 3, 5 = 0.3, 2.9). For monoatomic, diatomic, and Fibonacci chains, red dots 
represent 7 = 0.1, and blue dots represent 7 = 1. 

One can see that the occupation density of periodic chains is constant except 
at the edge, and the particle current reaches non-zero constanlRI for large n, and 



the constant depends on the parity of system size n for diatomic chains 
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the transport is ballistic. Red dots (7 = 0.1) show a small deviation since the 
interaction between the system and the reservoirs is relatively small. The Fibonacci 
chain shows large fluctuations in occupation and system size dependence of the 
particle current, though it is very robust against the change of 7. For small disordered 
chains (t = 3, 5 = 0.3), the particle current decreases linearly as a function of system 
size, and the occupation shows a linear profile. For disordered chains (t = 3, 5 = 
2.9), occupation profile has a kink shape and the particle current decreases nearly 
exponentially. 

For sufficiently small thermal and chemical gradients, the particle and the heat 
current, defined as Jq = Je — ^Jp (fi = (/Ul+/2r)/2), depend linearly on the external 
gradients asP®'^ 

J Q = L QQ A(3 - L QP (3An , 
J P = L PQ A(3 - LppPAfi , 

where Af3 = Pr — Pl and An = fiR — fip. The second law of thermodynamics imposes 
definite-positiveness of the matrix of Onsager coefficients L, which implies Lqq > 
and Lpp > 0, and if the dynamics is time-reversible, the Onsager 's reciprocity 
relation Lpq = Lqp holds. 

In Fig. [5] we consider the Fibonacci chains with = 3 and show the dependence 
of various properties of L on the other hopping parameter tp. Fig. |H(a) shows the 
tp dependence of all Onsager coefficients. One sees that all coefficients are positive, 
where we remark that possibility of negative off-diagonal elements were discussed 
inP 

Fig. fflb) shows the thermoelectric figure-of- merit ZT = LpgLgp/detlP for 
diatomic and Fibonacci chains. One sees that diatomic chains have larger ZT than 
the Fibonacci chains for the most parameter regimes. 

Fig. ID^c) shows the 7 dependence of \Lpq/Lqp — 1| for diatomic and Fibonacci 
chains. We see that the Onsager reciprocity is roughly linearly broken by increasing 
7 for the diatomic chains. It is very similar for the Fibonacci chains, however here 
the asymmetry has a cusp shape near 7 ~ 0.03, and one should be more careful with 
the choice of 7 in order to have (approximate) time-reversal symmetry. 

§5. Conclusions 

We have established a conservation law for particle number, and a conservation 
law for the sum of energy current along the chain and dissipation at the boudaries 
of the chain. The conservation laws are valid for generic one-dimensional chains 
connected to meso-reserovirs, which desribed by bilinear hamiltonian. 

As an application, we have studied the transport properties of monoatomic, 
diatomic, quasi-periodic, and disordered chains. We have observed wide fluctuations 
in the occupation and the system size dependence of the particle current for the 
Fibonacci chain. For the disordered chain, we have observed linearly or exponentially 
decreasing currents. The occupation was shown to have either linear profile or kink 
shape, respectively. 
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Fig. 1. (Color online) Panel (a) and (c) show position (j) dependence of occupation density rij 
for (a) monoatomic chain, and (c) diatomic chain. Panels (b) and (d) show n dependence of 
the particle current J p for (b) monoatomic chain, and (d) diatomic chain (red: 7 = 0.1, blue: 

7 = 1) 
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Fig. 2. (Color online) Panel (a) shows j dependence of occupation density rij for the Fibonacci 
chain. Panel (b) shows chain size (n) dependence of the particle current J p for the Fibonacci 
chain (red: 7 — 0.1, blue: 7 = 1) 
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Fig. 3. (Color online) Panel (a) and (c) show j dependence of occupation density rij for disordered 
chains. Panel (b) and (d) show n dependence of the particle current J p for disordered chains, 
(a) and (c) is for t — 3, S = 0.3, and (b) and (d) is for t — 3, S = 2.9 
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Fig. 4. (Color online) Panel (a) shows tp dependence of the Onsager coefficients (color code: Lqp 
- red, Lpp - black, Lqq - green, Lpq - blue) for the Fibonacci chains. Panel (b) shows tp 
dependence of ZT for diatomic (blue) and Fibonacci chains (red). Hoppings tA are chosen as 
tA = 3. For diatomic chains, they have exponentially small current outside if, £ [1, 7] indicated 
in the figure. Panel (c) shows the 7 dependence of asymmetry Lqp /Lpq — 1 for diatomic 
(blue) and Fibonacci (red) chains (dashed line indicates linear growth). tA and tp are chosen 
as tA = 3, tp — 6. 



